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We construct a nonlinear a model to describe a system of non-interacting electrons propagating 
in the presence of random magnetic flux. We find a term describing the long ranged logarithmic 
interaction between the topological density of the non-linear sigma model, and argue that this could 
give rise to a Kosterlitz-Thouless transition from the localized phase to a phase with power law 
correlations and continuously varying conductances. We provide a physical interpretation of our 
results in terms of the scattering of edge states of the magnetic domains in different regions. 
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The localization problem in two dimensions has been 
investigated extensively in the past twenty years . Re- 
cently there has been considerable interest in studying 
the problem of noninteracting two-dimensional electrons 
propagating in the presence of random magnetic flux. 
Physically this problem arises in the study of the quan- 
tum Hall systems near even denominator filling fractions 
and could also be related to holon propagation in 
chiral spin liquid Q and the gauge theory describing hole 
motion in an antiferromagnet [|6| . Single particle proper- 
ties, such as the one-body Green's function and the den- 
sity of states, have been addressed [|7|-[Tot. A more diffi- 
cult question to assess is that of localization. Naively, one 
would assume that the conventional unitary class sigma 
model description holds, in which case the beta function 
vanishes to the first order, but is non- vanishing and nega- 
tive to the second order Jl^ , indicating weak localization. 
For the random flux tight binding model, 
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with the phase e*"*' on each link uniformly distributed 
on the unit circle, numerical calculations in squares and 
long strip geometries are consistent with the fact that 
states are localized near the bottom of the band ||l^,|l3|,Q . 
However, as the energy approaches the band center, the 
localization length increases rapidly, and the interpre- 
tation of numerical data becomes problematic. In refs. 
P,p"3|P] it is was argued that states near the band center 
are extended, whereas ref. jl^ argues that all states are 
localized. In view of the rapid increase of the localization 
length and the limit of sample sizes in numerical calcu- 
lations, it is important to develop analytical methods to 
address the issue of localization. 

In this work, we construct a field theory to describe the 



random flux model. We envision a system with both ran- 
dom magnetic flux in addition to a random scalar poten- 
tial. Our starting point is the non-linear sigma model of 
the unitary ensemble which is appropriate for the present 
case with broken time reversal invariance. In the lattice 
model defined by equation (1), the random flux is cor- 
related only over the range of an elementary plaquette. 
However, for the field theory formulation it is more conve- 
nient to consider the opposite limit where the correlation 
length of the random magnetic field is large compared to 
the elastic mean free path. While these two models ad- 
dresse two different limits, we hope that the basic physics 
continuously interpolates between them. The basic de- 
grees of freedom of the field theory model are the matrices 
Q defined on the coset space of U{2m)/U{m) x U{m); in 
the localization problem we are interested in the limit of 
vanishing m. One remarkable feature of this field the- 
ory is the existence of the topological excitation called 
instantons. It was argued in ref. that in the presence 
of a uniform magnetic field, there is a term describing the 
coupling of the Hall conductance to the instanton density. 
This topological term is responsible for the existence of 
extended states necessary for the integer quantum Hall 
effect. In the case of random magnetic flux with no net 
average, such a topological term can not exist. However, 
it is perfectly sensible to think of a "local topological 
angle" associated with a given magnetic domain, which 
is coupled to the local instanton density. We are able 
to demonstrate that there is a finite stiffness to changes 
this local topological angle. This leads to an effective 
long ranged logarithmic interaction between the instan- 
ton density. We further argue that this long ranged in- 
teraction could lead to a Kosterlitz-Thouless transition 
between the free and confined phase of the instantons. 
In the former phase all states are localized whereas ex- 
tended states exist in the later case and there is gapless 
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excitations in the system which gives power law correla- 
tions. 

Physically, one can think of these gapless excitations 
as the edge states of the fluctuating magnetic domains. 
These extended edge states can exist over a distance 
larger than the localization length, just as in the case 
of the quantum Hall effect. These edge states are con- 
fined in the one dimensional network where the magnetic 
field changes sign. Since the magnetic field is randomly 
distributed with a zero average, this network always per- 
colate the entire system. Of extreme importance is the 
scattering of the edge states with different sense of rota- 
tion, or different chirality. In the phase where this scat- 
tering is relevant, a mass gap opens up and all states 
are localized. In the phase where the scattering is irrele- 
vant, the edge states remain gapless and carry currents. 
There are power law correlations and a finite, continu- 
ously varying conductance. 

We begin with the non-linear sigma model derived by 
Levine, Libby and Pruisken for the case of electrons 
propagating in a random scalar potential and an uniform 
magnetic field. The Lagrangian takes the following form: 
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C = -a,, Tr {d^^QY - -a^y Tr e^'Q d^Q d^Q (2) 

where Q is a matrix defined on the coset space of 
U{2m)/U{m) x U{m), and a^x and axy are the longi- 
tudinal and the Hall conductances respectively, defined 
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'^ly = {^ly) = - (Tr e^^H^G+H^G ))i„p 

< = {S^'xy) = *(Tr e^.r^H,(G+ - G-))i,„p (3) 

In the above formula, denotes the position operator 
and H^ denotes the gauge invariant momentum operator. 
G"*" and G~ are the advance and retarded one particle 
Green's functions respectively. (. . .)imp denotes impu- 
rity average whereas (. . .)imp denotes both quantum me- 
chanical trace and impurity average. The second term in 
the Lagrangian is called a topological term since for any 
smooth field configurations. 
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where n is an integer, called the Pontrjagin index. Dif- 
ferent n label different topological sectors with a given 
net instanton number. The presence of the topological 
term violates both time reversal and parity, and gives a 
handness to the two dimensional plane. Such a term is 
not allowed on symmetry grounds for the random flux 
problem with no net averaged magnetic field. 

cr^y is a one body operator that measures the magne- 
tization at the fermi energy, it receives only contribution 
from the edge states 15 . In the random field problem. 



its average vanishes, but its correlation due to the fluc- 
tuating magnetic domains is long ranged. To see this, we 
notice that its contribution to the topological term arises 
in a cumulant expansion of 

(expi jcPra^}^{r)p{r))=l + \ fdM»xy{r))p{r) + 
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d'rd\'{a'J^{ryj^{r'))p{r)p{r') + ... (5) 



where p{r) = Tr e^'^ Q d^Q d^Q is the local Pontrjagin 
density. In the presence of an uniform magnetic held, 
the second term in nonzero, reexponentialting it gives 
the topological term in the Lagrangian of equation (||). 
In this case, one has to reorganize the cumulant expan- 
sion into an expansion in a]^y — {cr^xy)- I* <^8,n be eas- 
ily shown that the correlation function of this operator 
is short ranged, therefore keeping higher order terms in 
the cumulant expansion would only give rise to irrelevant 
higher derivative terms. In the case of random flux with 
no net average, the second term in equation (||) vanish 
identically. One is left to calculate the two point correla- 
tion function of local a operators. This correlation func- 
tion can be calculated using a Ward identity relating it to 
the current current correlation function. Define the local 
current operator to be Jf_i{r) — iHp(r)(G"'" — G~). For 
the purpose of calculating DC transport properties one 
can assume that the current is transverse, i.e df^Jf^ = 0. 
Furthermore let us assume that the current operator can 
be expanded in terms of its moments. Under these two 
assumptions, 
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This is nothing but the familiar identity J = x M 
in electromagnetism, expressing current in terms of the 
magnetic moment, which in our case is just al!y. Since 
this is an operator identity, it can be used to derive the 
following identity between the correlation functions: 

{Mq)U-q)) = l{S,.uq'-qf.q.){iy'Iy{q)<y'Iy{-q))+o{q^) 

(7) 

In the limit of g — > 0, the transverse current current 
correlation of the left hand side approaches (Txx, therefore 
we obtain 



{^xyiq)<^xy{-q)) 



4:<7x 



(8) 



Since al.y is already a two body operator, its correlation 
function enters higher order cumulant expansion which 
will not be considered here. Inserting (|^) into equation 
(^) and reexponentialte the third term, one obtains the 
following expression for the field theory action of the ran- 
dom flux model: 
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Thus we see that while there is no topological term in 
the action for the random flux problem, there is a term 
describing the long ranged interaction of the topological 
density. 

The physical meaning of this nonlocal interaction term 
becomes transparent if one introduces an additional bo- 
son field to decouple the nonlocal term. The action be- 
comes local after this decoupling: 

£ = ^ Tr {d^Qf + — {d^cl^f + cj^ir) p{r) (10) 

The above equation is the central result of this paper, 
we propose this action to be the effective field theory 
model describing the random flux problem. The last term 
in the above equation is reminiscent of the topological 
term, except that the constant topological angle is now a 
space dependent scalar field, and there is a gradient term 
for the scalar field. The gradient term keeps (f> locally 
constant. The locally constant (j> field can be viewed as 
the topological angle in a given magnetic domain in the 
random flux problem. Another physical interpretation 
can be obtained when we use the fermionic representation 
of the 4> field. According to the bosonization rules, one 
scalar field can be represented by a Dirac fermi field, with 
the gradient term (5^(/)(r))^ replaced by the Dirac action 
"07^9^ "0- The equation of motion for the bose field 4> is 

dJ—d^A^P (11) 

The fermionized version of this equation is 

dJ^4'Tori)^P (12) 

This is the famous chiral anomaly equation that arises 
in the resolution of the U{1) problem in the quantum 
chromodynamics jl^. It is also similar to a model of 
holes moving in antiferromagnet background in 1-1-1 di- 
mensions . The spin degrees of freedom is described 
by the 0(3) nonlinear sigma model whereas the mass- 
less fermions describes the holes in the antiferromagnetic 
background. These two degrees of freedom are coupled 
exactly through the anomaly equation (|l2|) . At the classi- 
cal level, the two chiral components of the Dirac fermion 
are decoupled since the mass term is absent. However, 
this symmetry is broken at the quantum level by the 
presence of instantons. In the topologically trivial sec- 
tor, the integrated Pontrjagin density vanishes and the 
chiral charge is conserved. The presence of instantons 
violates the conservation of the chiral charge. 

Physically, one can view the two chiral components of 
the Dirac fermion as the gapless edge states of the fluc- 
tuating magnetic domains JTst . The right chiral fermion 



can be associated with the magnetic domain in the +z di- 
rection, and the left chiral fermion can be associated with 
the magnetic domain in the opposite direction. One can 
visualize the world line of the two chiral fermions as the 
current loops of the edge states. The edge states can be 
extended over a distance far larger than the localization 
length, thus they appear to be gapless at the classical 
level. However, there is tunneling between the two chiral 
components, where the tunneling event in our formula- 
tion is associated with the instanton, and the tunneling 
rate given by the instanton action. A crucial question is 
whether the chiral anomaly ( pi] ) and (|lj) gives an effec- 
tive mass for the Dirac fermion or if the fermions remain 
massless even in the presence of instantons. Formulated 
in another way, the question is whether the scattering 
between the two different types of edge states is relevant 
or irrelevant. 

One way to answer this question is to integrate out the 
Q fields in the nonlinear sigma model (^^ and obtain 
an effective theory for the boson field cj). One notices 
that there is a symmetry in the problem: if one shifts 
the 4> field by a constant 9, the gradient term remain 
unchanged. One picks up a contribution 

9 j(frp{r) = 16T:i9n (13) 

from the coupling of the cj) field to the topological den- 
sity. The above equality holds because of the quanti- 
zation of the topological charge (||). For 9 = ^, this 
shift contributes a factor of e'^'^" = 1 to the action, thus 
(j){r) — > (j){r) -|- i is an exact symmetry of the problem 
that must be respected by the effective Hamiltonian in 
(p. A mass term for the field, rn^cj)^ for example, is 
forbidden by this symmetry. Because of this symmetry, 
the effective action for (p can be expanded in a Fourier 
series, 

= — (5^0)' - 2y cos(167r0) + (14) 

^xx 

The higher harmonics term have higher scaling dimen- 
sions compared to the lowest harmonic term, and they 
can be neglected in the analysis of critical properties. 
Thus, general symmetry considerations severely restrict 
the low energy effective action for the field to be the 
sine-Gordon model. At this level, y is an unknown pa- 
rameter. While this effective action for the field is 
obtained from general symmetry considerations, it can 
also be derived explicitly from a dilute instanton gas ap- 
proximation. It is known that the nonlinear sigma model 
for the Q matrices supports instanton solutions that sur- 
vives the replica limit m — > fl^ . They make a contri- 
bution of y = Z? exp(— 47rcra;a;) to the action, where the 
exponent is the classical action of an instanton and D is 
the fluctuation determinant around the classical instan- 
ton. Assuming that these instantons can be treated as 
independent variables, one obtains the following for the 
effective action of the (f> field 
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the sums and integrals can be done explicitly to give the 
effective action ([l4|). 

The effective action ( p^ is the sine-Gordon model for 
scalar field, and has a well known Kosterlitz-Thouless 
phase transition. After a rescaling of the field ((> one ob- 
tains 
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this is identical to the effective Lagrangian for the XY 
model provided one makes the identification of Uxx 
with the inverse temperature 13—1/ kT and y with the fu- 
gacity of the XY vortices. The Kosterlitz-Thouless tran- 
sition occurs at = 2/ it. The conductance Uxx which 
enters the effective field theory is the "bare" conductance 
calculated within the self-consistent Born approximation. 
CT^j. divides the regime where the cosine term is relevant 
from the one in which the cosine term is irrelevant. For 
(^xx < CT^a;, one is in the high temperature phase where in- 
stantons form a free plasma and the system is disordered 
with a finite mass gap. However, for axx > tr^^j one is in 
the low temperature phase in which the instantons form 
tightly bound pairs. There is line of fixed points with 
power law correlations. The massive phase is identified 
with the localized phase in which the scattering between 
the two different chiral edge states are relevant, whereas 
the massless phase is identified with the phase where ex- 
tended states exists and the scattering between the edge 
states is irrelevant. The current is carried by these gap- 
less edge states and they are responsible for the finite 
continuously varying conductance in this regime. This 
picture is perfectly consistent with the recent numerical 
results of ref . |^,|l^,|| , where it was suggested that a line 
of fixed points might exist in a finite region around the 
center of the band. The conductance was found to be 
finite and continuously varying in this energy range. 

However, a word of caution concerning the dilute in- 
stanton gas approximation is in order here. In a model 
where the massless fermions are cougled to the gauge 
fields of the abelian Higgs model the dilute in- 

stanton gas approximation is valid because the instan- 
ton has a well defined size. Consequently the existence 
of the Kosterlitz-Thouless phase transition is on a firm 
ground. For the case of massless fermions coupled to non- 
linear sigma model, the instantons have arbitrary size 
and therefore they can overlap with each other which 
makes dilute instanton approximation hard to justify. 
This question was studied extensively in the context of 
the CP""i model [pl|-p3l. In the large n limit, Witten 



PH has shown that the instanton disappear after quan- 
tum fluctuations are taken into account, and the field 
is always massive. However, it is now well established 
that this large n result does not hold for n < 2. For the 
case of n = 2, the CP"~^ model is equivalent to the 0(3) 
nonlinear sigma model. From the analogy with quantum 
spins chains, it is known that this model has a second or- 
der phase transition when the topological angle is equal 
to TT, whereas the large n theory predicts a first order 
phase transition |2|]. For the case of n = 1, the CP'^~^ 
model is equivalent to the abelian Higgs model, where 
the instanton calculus is certainly correct. Based on this 
general arguments, Affleck ||2^ has argued that the large 
n result does not hold for n <2. It is thus plausible that 
the qualitative feature obtained by the instanton calculus 
is valid in this regime, and in particular, in the replica 
limit of n ^ that we are interested here. 

In conclusion we have constructed an effective field the- 
ory of the localization problem in random magnetic flux. 
This field theory has the same degrees of freedom as the 
nonlinear sigma model of unitary ensemble, but there 
is a long ranged logarithmic interaction of the topolog- 
ical density. This model can be reformulated in terms 
of a scalar field coupling to the topological density or 
equivalently, a Dirac fermi field obeying a chiral anomaly 
equation. These additional degrees of freedom can be 
physically interpreted as the edge states of the different 
magnetic domains. Our field theory model describes the 
tunneling between the edge states of different chirality 
in terms of the instantons of the nonlinear sigma model. 
The resulting effective field theory for the scalar field is 
shown to be a sine Gordon model from general symme- 
try considerations and we argue that such a model has a 
Kosterlitz-Thouless phase transition from localized states 
to extended states with power law correlations. This pic- 
ture is consistent with the recent numerical calculations 
of the random flux model in which a divergent localiza- 
tion length was found and a finite region of continuously 
varying conductance was identified. Our results have 
important consequences for the localizations physics by 
demonstrating the existence of a finite region of extended 
states in two dimensions, and deepen our understanding 
of the "Hall metal" phase around the even denominator 
filling fractions in the quantum Hall systems. 

The authors would like to acknowledge interesting and 
stimulating discussions with J. Chalker, D.H. Lee, H. 
Levine, S. Libby, N. Read, R. Shankar and A. Tikofsky. 
Part of this work was carried out at the Aspen center for 
physics during the workshop on the quantum Hall effect. 



[1] P. A. Lee and T. V. Ramakrishnan, Rev. Mod. Phys. 57, 
287, (1985). 

[2] V. Kalmeyer and S. C. Zhang, Phys. Rev. B46, 9889 



4 



(1992). 

[3] B. I. Halpcrin, P. A. Lcc and N. Read, Phys. Rev. B47, 
7312 (1993). 

[4] V. Kalmeyer D. Wei, D. Arovas, and S. C. Zhang, "Two- 
Dimensional Localization in the Presence of Random 

Flux and the Quantum Hall System at Even Denomi- 
nator Filling Fractions" , Phys. Rev. B (in press, October 
1993). 

[5] Kalmeyer and R. B. Laughhn, Phys. Rev. Lett. 59, 2095 
(1987); X. G. Wen, F. Wilczek and A. Zee, Phys. Rev. 
B39, 11413 (1989). 

[6] N. Nagaosa and P. A. Lee, Phys. Rev. Lett. 64, 2450 
(1990). 

[7] C. Pryor and A. Zee, Phys. Rev. B46, 3116 (1992). 

[8] B. L. Altshuler and L. B. loffe, Phys. Rev. Lett. 69, 2979 

(1992) . 

[9] G. Gavazzi, J. M. Wheatley, and A. J. Schofield, Phys. 

Rev. B47, 15170 (1993). 
[10] D. V. Khvoschenko and S. V. Meshkov, Phys. Rev. B47, 

12051 (1993). 
[11] S. Hikami, Prog. Theor. Phys. 107, 213 (1992). 
[12] T. Sugiyama and N. Nagaosa, Phys. Rev. Lett. 70, 1980 

(1993) . 

[13] Y. Avishai, Y. Hatsugai, and M. Kohmoto, Phys. Rev. 
B47, 9561 (1993). 

[14] H. Levinc, S. B. Libby, and A. M. M. Pruisken, Phys. 
Rev. Lett. 51, 1915 (1983); JVucJ. Phys. B240, 30 (1984); 
Nucl. Phys. B240, 49 (1984); JVucJ. Phys. B240, 72 
(1984); A. M. M. Pruisken in The Quantum Hall Effect, 
R. E. Prange and S. M. Girvin, eds. ( Springer- Verlag, 
1990). 

[15] P. Stfeda, J. Phys. C 15, L717 (1982). 

[16] For a review, see S. Coleman, Use of instantons, in As- 
pects of Symmetry (Cambridge, 1988). 

[17] R. Shankar, JVucJ. Phys. B330, 433 (1990). 

[18] This is similar with the identification of the chiral 
fermions with the edge states made in the recent liter- 
ature. Physically we are dealing the chiral fermions in 
two dimensional space, whereas the following literature 
studies the dynamics in 1 -|- 1 dimensional space time. X. 
G. Wen, Phys. Rev. Lett. 64, 2206 (1990); D.H. Lee and 
X.G. Wen, Phys. Rev. Lett. 66, 1765 (1991); M. Stone, 
Ann. Phys. 207, 38 (1991). 

[19] J. Kogut, Rev. Mod. Phys. 51, 659 (1979). 

[20] C. G. Callan Jr., R. Dashen, and D. J. Gross, Phys. Rev. 
D16, 2526 (1977). 

[21] E. Witten, JVucJ. Phys. B149, 285 (1979). 

[22] L Affleck, JVucJ. Phys. B305, 582 (1988). 

[23] L Affleck Phys. Rev. Lett. 66, 2429 (1991). 



5 



